Journal  of  Power  Sources  196  (2011)  10328-10335 


ELSEVIER 


Contents  lists  available  at  SciVerse  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


in 

SllbudtS 


Evaluation  of  commercial  lithium-ion  cells  based  on  composite  positive  electrode 
for  plug-in  hybrid  electric  vehicle  applications.  Part  I:  Initial  characterizations 

Matthieu  Dubarry3,  Cyril  Truchot3,  Mikael  Cugnet3,  Bor  Yann  yaw3’*,  Kevin  Geringb, 

Sergiy  Sazhinb,  David  Jamison5,  Christopher  Michelbacher5 

a  Hawaii  Natural  Energy  Institute,  SOEST,  University  of  Hawaii  at  Manoa  1 680  East-West  Road,  POST  1 09,  Honolulu,  HI  96822,  USA 
b  Energy  Storage  &  Transportation  Systems,  Idaho  National  Laboratory,  P.O.  Box  1625,  Idaho  Falls,  ID  83415-2209,  USA 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  20  June  2011 

Received  in  revised  form  18  August  2011 

Accepted  18  August  2011 

Available  online  25  August  201 1 


Keywords: 
Lithium-ion  battery 
Composite  electrode 
PHEV 

Cell-to-cell  variations 
Incremental  capacity 


Evaluating  commercial  Li-ion  batteries  presents  some  unique  benefits.  One  of  them  is  to  use  cells  made 
from  established  fabrication  process  and  form  factor,  such  as  those  offered  by  the  1 8650  cylindrical  config¬ 
uration,  to  provide  a  common  platform  to  investigate  and  understand  performance  deficiency  and  aging 
mechanism  of  target  chemistry.  Such  an  approach  shall  afford  us  to  derive  relevant  information  without 
influence  from  processing  or  form  factor  variability  that  may  skew  our  understanding  on  cell-level  issues. 
A  series  of  1.9  Ah  18650  lithium  ion  cells  developed  by  a  commercial  source  using  a  composite  positive 
electrode  comprising  {LiMn1/3Ni1/3Coi/302  +  LiMn204}  is  being  used  as  a  platform  for  the  investigation  of 
certain  key  issues,  particularly  path-dependent  aging  and  degradation  in  future  plug-in  hybrid  electric 
vehicle  (PHEV)  applications,  under  the  US  Department  of  Energy’s  Applied  Battery  Research  (ABR)  pro¬ 
gram.  Here  we  report  in  Part  I  the  initial  characterizations  of  the  cell  performance  and  Part  II  some  aspects 
of  cell  degradation  in  2C  cycle  aging.  The  initial  characterizations,  including  cell-to-cell  variability,  are 
essential  for  life  cycle  performance  characterization  in  the  second  part  of  the  report  when  cell-aging  phe¬ 
nomena  are  discussed.  Due  to  the  composite  nature  of  the  positive  electrode,  the  features  (or  signature) 
derived  from  the  incremental  capacity  (IC)  of  the  cell  appear  rather  complex.  In  this  work,  the  method 
to  index  the  observed  IC  peaks  is  discussed.  Being  able  to  index  the  IC  signature  in  details  is  critical  for 
analyzing  and  identifying  degradation  mechanism  later  in  the  cycle  aging  study. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Rechargeable  lithium-ion  batteries  (LIB)  are  viable  choices  for 
numerous  energy  storage  and  mobile  power  applications.  Most 
notably  in  the  transportation  sectors,  plug-in  hybrid  electric  vehi¬ 
cles  (PHEVs)  are  of  great  interest  in  using  LIB  for  propulsion,  as 
promised  by  its  prospect  in  high  energy  efficiency  and  potential 
for  long  driving  range.  However,  reliable  deployment  of  LIB  in 
powertrains  remains  very  challenging,  and  more  understanding  on 
battery  performance  and  aging  mechanisms  including  cell  degra¬ 
dation  under  conditions  specific  to  PHEV  applications  to  achieve 
greater  efficiency  is  desired.  For  a  battery  to  reach  that  goal,  both 
high  energy  and  high  power  features  are  required. 

Extensive  testing  of  targeted  cell  chemistries  for  PHEV  appli¬ 
cations  is  necessary  to  understand  each  new  generation  of  battery 
chemistries  and  their  performance.  Several  generations  of  LIBs  have 
been  tested  in  the  DOE  national  laboratories  [  1  -3  ]  using  specifically 
designed  protocols  [4-7].  In  this  work,  the  targeted  cell  chemistry 
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uses  a  positive  electrode  (PE)  active  material  (PEAM)  comprising 
{LiMn1/3Ni1/3Co1/302  +  LiMn204}  composite  and  a  negative  elec¬ 
trode  (NE)  made  of  graphite.  The  concept  of  composite  electrodes 
is  attracting  attention  in  recent  years  because  it  promises  the 
combination  of  merits  of  several  active  materials  into  a  hybrid 
electrode  to  allow  performance  optimization  [8,9].  For  instance, 
LiMn204  (spinel)  is  a  promising  PEAM  with  a  high  electrode  poten¬ 
tial  against  the  Li  metal  reference,  excellent  rate  capability,  and 
relatively  low  cost  for  production;  however,  it  exhibits  a  rela¬ 
tively  low  capacity,  ca.  140-1 50  mAh  g_1,  and  suffers  degradation 
from  manganese  (Mn2+)  dissolution  in  the  electrolyte  and  the 
subsequent  deposition  on  the  NE  in  the  charge  regime  to  cat¬ 
alyze  electrolyte  decomposition,  especially  at  high  temperatures 
[9-11].  LiMn1/3Ni1/3 Co^i302  is  another  promising  PEAM  with  a 
higher  capacity,  around  180mAhg-1,  but  worse  rate  capability 
than  the  spinel.  The  composite  composition  is  an  ultimate  com¬ 
promise,  which  could  provide  an  electrode  that  exhibits  optimized 
capacity,  rate  and  power  capability,  in  contrast  to  more  conven¬ 
tional  PEAMs,  such  as  alumina-coated  LiCo02  (ca.  170  mAh  g-1 
[12]),  LiNi02  (ca.  170-lSOmAhg-1  [13,14]),  LiNi*Coi_x02  (ca. 
150-160 mAh g-1  [15]),  and  LiNio.8Coo.15Alo.05O2  (ca.  ^OmAhg-1 
[16]).  As  an  additional  benefit,  it  has  been  found  that  the  addition 
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of  the  doped  cobalt  oxides  prevents  the  dissolution  of  man¬ 
ganese;  therefore,  improving  the  capacity  retention  of  LiMn204 
upon  cycling  [9-1 1,17-21  ]. 

Through  extensive  testing  effort,  HNEI  developed  several  useful 
electrochemical  techniques  in  recent  years  allowing  the  character¬ 
ization  of  cell-to-cell  variations  within  a  batch  [22,23]  as  well  as 
the  understanding  of  how  to  employ  such  in  situ  electrochemi¬ 
cal  techniques  [24-28]  to  improve  the  ability  of  extracting  detailed 
information  from  data  obtained  from  commonly  used  test  protocols 
in  battery  performance  evaluations  and  degradation  studies. 

In  the  initial  conditioning  and  characterization  tests  (ICCT),  a 
modified  version  of  test  protocols  and  procedures  was  derived  from 
those  defined  in  the  USABC  Electric  Vehicle  Battery  Test  Procedures 
Manual  [29].  The  intent  of  the  ICCT  is  to  establish  the  baseline 
performance  matrix  for  the  cell  chemistry,  including  cell-to-cell 
variations.  The  cell-to-cell  variations  are  critical  to  us  in  several 
aspects  [22,23],  including: 

(1 )  To  quantify  the  cell  quality  in  performance. 

(2)  To  quantify  such  variability  in  performance  in  multi-cell  appli¬ 
cations. 

(3)  To  establish  the  confidence  level  for  interpreting  test  results 
from  a  large  matrix  of  cells  that  will  undergo  various  test  con¬ 
ditions. 

(4)  To  estimate  possible  errors  in  model  accuracy. 

(5)  To  estimate  possible  escalation  of  errors  in  subsequent  aging 
tests. 

Another  objective  of  the  ICCT  is  to  establish  the  capability  for  cell 
diagnosis  with  regard  to  degradation;  particularly,  the  relationship 
of  rest  cell  voltage  (RCV)  vs.  state  of  charge  (SOC)  [24]  and  the  basic 
features  to  allow  incremental  capacity  analysis  (ICA). 

Subsequently  various  duty  cycle  and  thermal  aging  tests  were 
conducted  in  Idaho  National  Laboratory  (INL)  and  University  of 
Hawaii  (UH),  respectively.  The  second  part  of  this  report  is  to 
describe  how  to  develop  an  in-depth  understanding  of  the  cell 
aging  process  to  enable  identification  of  cell  degradation  mech¬ 
anism^),  so  we  can  eventually  quantify  the  contribution  from 
each  mechanism  in  the  aging  process.  This  understanding  will 
enable  the  establishment  of  a  quantitative  physical  model  describ¬ 
ing  the  degradation  process  for  the  targeted  chemistry  with  an 
ability  to  accommodate  any  variance  in  cell  degradation  due  to  path 
dependence  [30]. 

2.  Experimental 

Ten  1 .9  Ah  1 8650  commercial  lithium-ion  cells  were  used  in  the 
study.  The  cells  were  surveyed  by  weight  and  open  circuit  voltage 
(OCV)  measurements  as  received  and  then  subjected  to  the  ICCT 
using  a  Maccor  4300®  system.  All  cells  were  tested  with  laboratory- 
made  racks,  and  our  experiences  indicate  that  it  is  important  to 
pay  attention  to  the  cable  and  its  contact  with  the  cell  to  minimize 
variability  in  the  configuration  of  the  test  harnessing  in  order  to 
achieve  better  consistency  in  test  results  and  improve  fidelity. 

Upon  the  initial  survey,  the  cells  were  subjected  to  a  condition¬ 
ing  process  using  a  series  of  charge  and  C/2  discharge  regimes, 
which  was  terminated  when  the  capacity  was  stabilized  within 
±0.2%  between  two  consecutive  cycles.  The  cells  can  be  condi¬ 
tioned  within  3-6  cycles  typically.  Unless  specified,  all  cells  were 
recharged  using  a  constant  current  (CC)  at  C/2  step  to  4.2  V,  fol¬ 
lowed  by  a  constant  voltage  (CV)  step  at  4.2  V  until  the  termination 
current  of  C/25  was  reached.  Unless  specified,  the  discharge  cutoff 
voltage  is  generally  2.6  V.  A  4-h  rest  between  charge  and  discharge 
regimes  was  imposed  to  determine  the  RCV.  The  cells  were  then 
subjected  to  two  discharge  regimes  at  C/25  and  C/5,  after  being 


fully  recharged  by  the  CC-CV  steps  each  times  to  determine  the 
corresponding  capacity,  respectively.  The  C/n  rate  represents  the 
discharge  current  calculated  from  the  nominal  capacity  Qnom  of 
1 .9  Ah  divided  by  the  discharge  duration  n  in  hours  (h).  The  capacity 
at  a  specific  C/n  rate  is  denoted  as  Qn-  The  polarization  resistance 
of  the  cell  is  calculated  from  the  difference  in  IR  drops  among  C/25, 
C/5  and  C/2  using  Ohm’s  law;  i.e.  A  V=  RAI,  assuming  that  the  polar¬ 
ization  in  the  cell  follows  a  pseudo-Ohmic  relationship  below  C/2. 
After  these  tests,  the  cells  were  recharged  to  50%  SOC  for  storage  at 
-10  °C  in  a  freezer. 

A  “nominal  sample  cell”  or  NSC  was  selected  from  the  batch 
and  subjected  to  additional  cycling  at  C/25,  C/5,  C/2, 1C,  2C,  and  5C 
to  determine  its  performance  characteristics.  In  these  cycles,  the 
cell  was  discharged  and  charged  at  specific  rates  until  the  cutoff 
condition  was  reached.  The  cell  was  set  for  4-h  rest  to  measure 
the  RCV  at  EOC  and  EOD  for  each  given  rate.  Remnant  capacity 
was  measured  following  the  previous  regime  with  a  C/25  dis¬ 
charge  or  charge  step  to  the  respective  cutoff  condition,  followed 
by  another  4-h  rest.  The  RCV  measured  at  the  end  of  this  rest 
shall  provide  the  SOC  of  the  starting  point  for  the  next  regime, 
either  noted  as  the  beginning-of-charge  (BOC)  or  as  the  beginning- 
of-discharge  (BOD),  when  appropriate.  To  determine  the  SOC,  we 
used  a  method  described  in  our  previous  work  [24-28],  which 
uses  the  pseudo-OCV  vs.  SOC  (denoted  as  ps-OCV  =  f(SOC))  curve 
for  this  chemistry.  The  ps-OCV  =  f(SOC)  curve  is  derived  as  follows: 
(1)  normalize  the  C/25  charge  and  discharge  (Vcell  =  f(Q))  curves  to 
VCeii  =  f(Q/Q25)  curves,  respectively,  for  the  NSC.  (2)  Average  the  Vcell 
of  the  Vcell  =  f(Q/Q25)  curves  in  the  charge  and  discharge  regimes. 
The  ps-OCV  =  f(SOC)  curve  has  been  verified  to  be  universal  within 
this  batch.  The  SOC  at  any  given  condition  can  be  determined  by 
inferring  the  corresponding  RCV  on  the  ps-OCV  =  f(SOC)  curve. 

The  incremental  capacity  (-dQ/dV=f(Vcell))  curves  are  derived 
from  correlating  the  amount  of  incremental  Q.  (in  Ah)  associated 
with  a  successive  voltage  step,  typically  in  4  mV,  in  the  charge  or 
discharge  curve  to  Vcell.  The  fidelity  of  the  -dQ/dV=f(Vcell)  curve 
is  sensitive  to  the  resolution  of  the  data  acquisition  system.  We 
typically  record  data  in  2  mV  intervals. 

3.  Results 

3.1.  Cell-to-cell  variations 

Fig.  1  displays  (a)  the  weight  and  (b)  initial  OCV  of  the  ten  cells 
as  received.  According  to  the  manufacturer’s  specification,  the  cell 
should  weigh  about  43  g.  All  the  cells  are  lighter  than  the  nomi¬ 
nal  weight  with  an  average  at  42.725  ±  0.058  g  (±0.13%).  The  small 
variations  on  weight  suggest  that  the  manufacturing  process  is 
quite  consistent.  The  average  initial  OCV  recorded  is  3.709  ±  0.001  V 
(±0.0(2)%),  which  corresponds  to  38.2  ±  0.2%  SOC  according  to  the 
ps-OCV  =  f( SOC)  curve  [24-28].  Thus,  all  ten  cells  have  almost  the 
same  SOC  as  received. 

Fig.  2(a)  shows  the  capacity  variations  obtained  from  C/25  and 
C/2  discharge  regimes.  The  capacity  measured  from  each  of  the 
three  C/n  rates  (C/5  is  not  shown)  is  quite  consistent  among  the 
cells,  further  demonstrating  the  consistency  of  the  cell  quality  in 
performance.  At  C/2  the  cells  retain  1.888  ±0.003  Ah  (±0.2%)  on 
average,  whereas  at  C/25  a  little  more  capacity  is  available  with  an 
average  of  1 .949  ±  0.003  Ah  (±0.2%).  The  cells  are  rated  at  1 .9  Ah  at 
C/2;  therefore,  all  the  cells  are  slightly  short  of  the  rated  capacity. 
It  is  also  worth  noting  that  for  rates  equal  to  or  lower  than  C/5  all 
cells  are  able  to  provide  the  rated  capacity. 

Fig.  2(b)  presents  the  RCV  at  the  end-of-discharge  (EOD)  of  C/2 
and  C/25.  As  expected,  the  lower  the  rate,  the  lower  the  RCV  at 
the  EOD.  At  the  EOD  of  C/25,  the  cell  relaxed  to  2.975  ±  0.005  V 
(±0.2%)  on  average,  whereas  at  C/2  it  is  3.364  ±  0.009  V  (±0.3%). 
The  consistency  in  the  RCV  at  the  end-of-charge  (EOC)  prior  to  the 
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Fig.  1.  (a)  Weight  and  (b)  as-received  OCV  for  the  ten  commercial  18650  cells. 
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Fig.  2.  (a)  Capacities  at  C/2  and  C/25  and  (b)  rest  cell  voltage  (RCV)  at  the  end-of- 
discharge  (EOD)  for  the  ten  cells. 

C/2  and  C/25  discharges  is  even  better,  as  the  deviation  is  about 
±0.002  V  (±0.1%)  on  average. 

Fig.  3  presents  the  polarization  resistance  derived  from  the  dis¬ 
charge  curves  for  the  ten  cells.  The  average  polarization  resistance 
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Fig.  3.  Distribution  of  polarization  resistance  in  the  ten  cells. 


Capacity  (Ah) 

Fig.  4.  Charge  and  discharge  curves  for  cell  #1  RPT  at  rates  ranging  from  C/25  to  5C. 


is  67.1  ±3.8m£2  (±5.7%),  which  contributes  to  little  disparities  in 
the  capacity  among  the  cells. 

3.2.  Initial  characterization  on  nominal  sample  cell 

Fig.  4  presents  the  Vcell  =  f(Qj  curves  for  various  charge  and  dis¬ 
charge  regimes  in  NSC  during  the  reference  performance  test  (RPT). 
It  should  be  noted  that  the  cell  has  been  subjected  to  remnant 
capacity  measurement  at  the  end  of  each  regime;  therefore,  the 
starting  SOC  is  the  same  for  all  rates  at  a  given  regime.  In  Fig.  5(a), 
the  capacity  obtained  from  various  charge  and  discharge  regimes 
are  reported.  In  the  discharge  regime,  the  capacities  for  rates  from 
C/5  to  2C  follow  the  Peukert  law  Q.=  /k  x  t;  where  k  is  the  Peukert 
coefficient  and  t  the  nominal  discharge  time  (in  h).  The  Peukert 
coefficient  k  =  1.019  was  determined,  which  indicates  that  the  cell 
has  reasonable  rate  capability.  Flowever,  the  Q25  appears  to  be 
lower  than  that  extrapolated  from  the  Peukert  equation  (1.949  Ah 
vs.  1.998  Ah),  suggesting  that  the  cell  maximum  capacity  in  this 
voltage  range  is  reached.  The  capacity  at  5C  is  below  the  value 
extrapolated  from  the  Peukert  equation.  The  Peukert  law  primarily 
takes  into  account  the  capacity  variability  originated  from  varia¬ 
tions  in  mass  transfer  in  the  porous  electrodes.  It  does  not  take 
into  account  any  polarization  effect.  It  assumes  that  all  electrode 
reactions  were  completed  at  a  given  rate.  Judging  from  the  voltage 
response  in  the  5C  discharge  regime  as  shown  in  Fig.  4,  it  appears 
that  the  cutoff  voltage  has  been  reached  prior  to  the  completion  of 
the  electrochemical  reactions  at  the  surface  of  the  electrodes.  We 
further  verified  that  the  cutoff  was  indeed  reached  prematurely 
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before  the  knee  that  is  typical  for  the  EOD.  This  premature  cutoff 
explains  why  the  5C  capacity  does  not  follow  the  Peukert  law. 

In  the  charge  regime,  due  to  the  same  constraint  and  prema¬ 
ture  cutoff,  the  rate-dependent  charge  return  did  not  follow  the 
Peukert  law  either.  It  is  understandable  that,  due  to  the  stability 
window  of  the  electrolyte  and  the  nominal  cell  voltage,  which  is 
close  to  3.7  V,  the  cell  has  a  1.0  V  margin  in  the  discharge  regime; 
yet,  only  0.5  V  in  the  charge  regime,  to  accommodate  polarization. 
It  is  therefore,  similar  to  most  of  the  high-voltage  chemistries,  that 
this  cell  chemistry  is  susceptible  to  overcharging  under  high-power 
pulse  charging  conditions  (such  as  in  regenerative  breaking). 

Conventionally,  the  power  capability  of  the  cell  should  be  deter¬ 
mined  by  constant  power  discharge  tests.  The  assessment  of  power 
capability  could  be  approximated  by  the  power  requirements 
derived  from  constant  current  tests  by  multiplying  the  average 
voltage  with  the  C/n  rate  current.  Such  an  estimate  of  power  capa¬ 
bility  often  results  in  an  underestimate  from  the  one  by  constant 
power  test  because  in  the  case  of  constant  current  discharge  the 
polarization  increases  more  aggressively  toward  the  EOD  and  leads 
to  a  lower  power  amplitude.  Nevertheless,  this  approximation 


provides  a  reasonable  estimate  of  the  trade-off  between  power  and 
energy  for  this  chemistry,  as  revealed  in  the  Ragone  plot  in  Fig.  5(b). 
The  cell  is  able  to  deliver  125  Whkg-1  at  650  W  kg-1,  indicating  a 
good  power  capability  of  this  chemistry.  In  the  charge  regime,  how¬ 
ever,  the  energy  absorbable  by  the  cell  drops  below  100 Whkg-1 
at  250  W  kg-1,  suggesting  that  the  chemistry  cannot  handle  high 
charge  rates  for  fast  charging  or  regenerative  braking  over  a  con¬ 
siderable  range  of  high  SOC. 

4.  Discussion 

4.1.  Cell-to-cell  variations 

We  recently  demonstrated  that,  using  common  testing  proce¬ 
dures  and  techniques  to  analyze  the  conditioning  cycles  of  100  cells 
provided  by  a  vendor,  three  attributes  were  derived  to  characterize 
the  cell  variations:  the  amount  of  active  material  in  the  cell,  polar¬ 
ization  resistance,  and  other  kinetic  factors  that  confine  the  rate 
capability  [22,23].  The  amount  of  active  material  dictates  the  max¬ 
imum  capacity  in  a  cell.  The  polarization  resistance  dictates  the  cell 
voltage  under  polarization,  thus  the  portion  of  capacity  and  power 
deliverable  under  a  certain  rate  and  cutoff  conditions.  Finally,  the 
rate  capability,  the  capability  of  a  cell  in  handling  discharge  rate 
variations,  is  dictated  by  the  additional  kinetic  factors,  as  expressed 
by  the  Peukert  coefficient. 

The  amount  of  active  material  in  the  cell  can  be  determined 
by  the  capacity  released  at  a  very  slow  rate  such  as  C/25  or  less; 
or,  from  a  process  that  determines  the  quantity  called  “capacity 
ration,”  in  Ah  %SOC-1,  which  is  derived  from  the  capacity  at  a 
given  rate  and  the  RCVs  measured  at  EOC  and  EOD  that  denote 
the  SOC  range  involved  in  delivering  the  said  capacity.  From  the 
RCVs,  the  SOC  range  can  be  determined  from  the  ps-OCV  =  f(SOC) 
curve  for  the  given  rate.  Dividing  the  capacity  with  this  SOC  range, 
the  “capacity  ration”  is  thus  calculated  [22,23].  It  is  important  to 
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Fig.  6.  (a)  Capacity  ration  and  (b)  Peukert  coefficient  variations  for  the  10  cells. 
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Fig.  7.  (a)  Evolution  of  the  RCV  with  rate,  (b)  ps-OCV  =  f(SOC)  curve  calculated  from  C/25  charge  discharge  curves,  and  (c)  the  correspondence  of  SOC  range  between  BOD  and 
EOD  and  capacity  at  various  rates. 


recognize  that  the  two  methods  depicts  the  same  concept  that  the 
truly  meaningful  SOC  can  only  be  determined  by  electrochemi¬ 
cal  equilibrium,  in  which  the  range  of  SOC  should  coincide  with 
the  capacity  at  any  rate;  and,  such  a  “capacity  ration”  associates 
with  and  defines  the  content  of  the  active  materials  in  the  cell. 
It  is  also  important  to  note  that  the  capacity  ration  is  thus  rate 
independent. 

The  additional  kinetic  constraints  other  than  the  polarization 
resistance  are  usually  introduced  in  the  manufacturing  pro¬ 
cesses,  e.g.  inhomogeneity  arisen  from  differences  in  the  electrode 
thickness,  additive  ratio,  and/or  electrode  porosity  [31 ,32].  We  pos¬ 
tulated  that  these  factors  would  influence  the  effective  utilization 
of  active  material  in  general  and  thus  induce  variations  in  the  rate 
capability  performance  of  a  given  cell,  which  can  be  represented  by 
variations  in  the  Peukert  coefficient.  The  variation  in  the  rate  capa¬ 
bility  can  be  estimated  from  the  capacities  measured  in  the  ICCT 
(thus,  at  C/25,  C/5  and  C/2). 

Fig.  6  presents  the  values  of  (a)  the  capacity  ration  and  (b) 
Peukert  coefficient  for  the  ten  cells.  The  average  capacity  ration 
is  19.66  mAh  %SOC-1  and  its  deviation  among  the  cells  is  ±0.2%. 
The  average  Peukert  coefficient  is  1 .01 3  ±  0.1  %.  The  small  deviation 
observed  indicated  that  the  intrinsic  quality  of  the  cell  assembly  is 
rather  consistent,  as  shown  by  the  rate  capability.  The  polarization 


resistance  (shown  in  Fig.  4)  varies  between  63  and  72  m^  (about 
6%  in  deviation). 

The  overall  performance  variations  among  the  cells  are  rel¬ 
atively  insignificant,  compared  to  other  results  of  different  cell 
chemistries  we  measured  so  far.  The  discharge  curves  among  the 
cells  can  be  superimposed  upon  one  another  even  at  C/2,  which 
indicates  that  the  cell  technology  may  support  intended  PFIEV  well 
in  pack  applications  and  that  the  results  of  the  NSC  can  be  extended 
to  others  in  the  batch. 

4.2.  Initial  characterization  on  nominal  sample  cell 

Fig.  7(a)  presents  the  evolution  of  RCV  for  each  charge  and  dis¬ 
charge  regime  as  a  function  of  C  rate.  Since  the  same  remnant 
capacity  measurement  and  a  4-h  rest  were  imposed  at  the  end  of 
each  regime,  RCVs  at  BOC  or  BOD  are  rather  independent  of  C  rate; 
thus,  all  regimes  start  with  almost  the  same  SOC  at  BOC  or  BOD, 
respectively.  However,  the  RCVs  at  EOC  or  EOD  vary.  The  EOC  RCVs 
follow  a  linear  relationship  with  C  rate,  suggesting  that  the  cut¬ 
off  was  primarily  determined  by  polarization.  The  EOD  RCVs  are  a 
strong  function  of  C  rates  below  1C  and  approach  a  steady  value 
after  2  C.  To  understand  this  behavior,  we  refer  the  RCVs  to  the  ps- 
OCV  =f(SO  C)  curve  as  shown  in  Fig.  7(b).  A  voltage  plateau  exists 


M.  Dubarry  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  10328-10335 


10333 


a 


b 


Voltage  (V) 


JO 

u. 

3 

O 


Positive  electrode  Potential  vs.  Li/Li+  Lithium  content 


© 

0 

(D 

© 

© 

© 

©  : 

:  ©:  © 

© 

© 

© 

© 

© 

© 

©■ 

*  : 

:  *  ! 

:  * 

* 

* 

* 

* 

* 

*  : 

©  ; 

!©;© 

© 

© 

© 

O 

O 

O 

©  ! 

© 

O 

© 

© 

O 

Negative  electrode  Lithium  content 


Fig.  8.  (a)  Incremental  capacity  (IC)  signature  of  cell  #1  at  C/25  and  (b)  schematic 
presentation  of  the  cell  reactions  in  the  IC  peak  index. 


around  3.5  V,  which  explains  why  the  RCVs  remain  rather  constant 
for  2C  and  5C,  whereas  the  capacity  varies  substantially.  The  pres¬ 
ence  of  the  plateau  on  the  ps-OCV  =  f(SOC)  curve  indicates  that  both 
electrodes  undergo  phase  transformations.  These  reactions  are  dis¬ 
cussed  in  details  in  Fig.  8  and  this  plateau  is  shown  to  be  associated 
with  a  graphite  staging. 

Fig.  7(c)  presents  the  comparison  between  the  capacities  nor¬ 
malized  with  Q25  in  the  test  and  calculated  from  the  RCVs  (the 
corresponding  SOCs  inferred  from  Fig.  7(b)).  The  excellent  agree¬ 
ment  supports  the  validity  of  our  approach  in  using  ps-OCV  =  f(SOC) 
curve  inference  to  determine  SOC.  Several  implications  arise  from 
such  validity.  The  mixed  oxide  cathode  and  carbonaceous  anode  in 
the  cell  chemistry  should  follow  the  core-shell  model  and  Li  diffu¬ 
sion  from  the  surface  to  the  core  of  the  grains  is  the  rate-limiting 
step  in  the  cell  reaction  and  relaxation  process.  The  rate  of  diffusion 
is  probably  comparable  to  C/8,  as  suggested  by  the  achievement  of 
the  maximum  capacity  at  this  rate  (as  estimated  from  the  Peukert 
equation).  Therefore,  at  rates  lower  than  C/8  all  the  active  materials 
are  involved  in  the  electrochemical  reactions  to  achieve  the  max¬ 
imum  capacity,  whereas  at  rates  higher  than  C/8  the  core  of  the 
grains  remains  deficit  in  Li  up-taking;  thus,  only  partial  capacity  is 
delivered. 

The  targeted  cell  chemistry  in  this  work  employs  a  compos¬ 
ite  PE  (c-PE)  comprising  {LiMn1/3Ni1/3Co1/302  +  LiMn204}  and  a 
graphite  NE.  Thus,  a  complex  cell  incremental  capacity  (IC)  signa¬ 
ture  is  expected.  Although  a  more  definitive  derivation  of  the  IC 
signature  should  be  performed  from  half-cell  experiments  for  both 
electrodes  in  reference  to  Li  metal,  it  is  unfortunately  impossible 
for  us  to  conduct  such  experiments.  Therefore,  to  infer  such  com¬ 
plex  electrode  reactions  in  reference  to  the  IC  signature  can  only  be 
synthesized  from  literature  reports.  If  the  two  active  materials  in 
the  c-PE  do  not  interfere  with  each  other,  we  infer  that  IC  signals 
of  each  active  material  in  the  c-PE  shall  be  convoluted  with  those 
of  the  graphite  NE,  resulting  in  a  series  of  distinct  IC  peaks.  Fig.  8(a) 
shows  the  IC  behavior  of  the  NSC  discharged  at  C/25  and  Fig.  8(b) 
displays  the  peak  index  based  on  electrode  reactions  inferred  from 
those  reported  in  the  literature  [32-36]. 

In  the  PE  potential  range  vs.  a  lithium  reference  electrode  in  this 
study,  LiMn1/3Ni1/3Co1/302  presents  a  steady  potential  decrease 
from  4.2  V  to  3.75  V  in  the  solid  solution  region  in  the  hexagonal 


phase.  Further  Li  intercalation  in  the  crystal  structure  leads  to  a 
phase  transformation  to  the  rock  salt  phase  and  exhibits  a  poten¬ 
tial  plateau  at  3.75  V  [35,36].  In  the  case  of  LiMn204  two  potential 
plateaus  at  4.1  V  and  3.95  V  are  reported  for  two  phase  transfor¬ 
mations  [34].  Accordingly,  two  IC  peaks,  a  sharper  one  at  4.1  V 
and  another  with  a  broader  asymmetric  shape  at  3.95  V;  are  thus 
expected.  Summarized  from  these  literature  reports,  we  can  then 
infer  the  IC  curve  for  the  c-PE  with  four  IC  peaks  as  follows:  the  first 
peak  (a  broad  band  from  4.2  V  to  3.75  V)  is  the  one  for  the  solid  solu¬ 
tion  region  in  LiMn1/3Ni1/3Co1/302  (CD ).  Within  this  broad  potential 
range,  the  two  peaks  (4.1  V  and  3.95  V)  in  LiMn204  (CD  and  (D )  shall 
intrude  before  the  phase  transformation  in  LiMn1/3Ni1/3Co1/302  (© 
)  takes  place  at  3.75  V. 

In  the  NE,  when  discharged,  the  lithiated  graphite  undergoes 
five  distinct  phase  transformations,  as  denoted  by  O  to  © ,  according 
to  the  literature  [32].  Fig.  8(b)  summarizes  the  reactions  in  the  two 
electrodes  on  a  unified  potential  scale.  In  order  to  index  the  IC  peaks 
by  their  distinctive  reactions  in  the  electrodes,  we  use  the  following 
notation  as  convention:  (1 )  an  IC  peak  is  a  unique  representation  of 
a  NE  reaction  convoluting  with  a  corresponding  one  on  the  positive, 
(2)  the  NE  reaction  is  denoted  by  a  black  (filled)  numbering,  while 
the  positive  with  an  open  (unfilled)  one,  (3)  using  “*”  to  denote 
the  convolution  between  the  negative  and  positive  IC  signals  that 
results  in  the  unique  IC  peak  in  the  cell.  Thus,  the  peak  indexed  as 
O  *(D  represents  reaction  O  in  the  NE  convoluting  with  reaction  (D 
in  the  PE. 

From  the  indexing,  we  found  that  some  features  of  the  IC  peaks 
can  be  used  as  status  indicators  for  an  electrode  (or  its  constituents), 
if  a  potential  plateau  is  present  on  the  opposite  electrode.  For 
instance,  peak  O  *(D  may  reveal  the  status  of  the  state  of  LiMn204, 
since  O  is  on  a  potential  plateau;  likewise,  peak  ©  *©  is  an  indicator 
for  LiMn!  /3  Nq  /3  /3  02 ,  and  peak  ©  *©  for  the  graphite  electrode.  In 
addition,  the  position  and  area  of  a  few  specific  IC  peaks  can  be  used 
to  estimate  the  composition  ratio  between  the  two  PE  constituents. 
By  examining  the  IC  peak  area  associated  with  © ,  which  accounts 
for  1  /  3  the  total  cell  capacity,  as  well  as  half  the  LiMn!  /3  N^  /3  Co!  /3  02 
capacity,  we  thus  conjecture  that  the  c-PE  composition  likely  com¬ 
prise  2/3  ofLiMn1/3Ni1/3Co1/302  and  1/3  of  LiMn204.  Similarly,  by 
examining  the  capacity  associated  with  peaks  O  on  the  NE  we  can 
estimate  the  loading  ratio  between  PE  and  NE  (after  the  initial  SEI 
formation  on  the  NE).  Ideally,  O  accounts  for  half  of  the  capacity 
of  a  typical  graphite  NE,  since  it  corresponds  to  the  transformation 
from  LiCi2  to  LiC6.  A  useful  approach  to  determine  the  electrode¬ 
loading  ratio  is  to  study  the  capacity  measured  at  low  rates,  in 
which  the  liberation  of  capacity  undergoes  with  a  minimal  influ¬ 
ence  from  kinetics  while  delivering  its  maximum  Ah  content  in  the 
active  materials  in  the  cell.  For  instance,  the  Q.  associated  with  O 
is  41.5%  of  the  cell  capacity.  Since  the  graphite  is  near  complete  Li 
depletion  at  the  EOD  (as  conjectured  by  the  presence  of  peak  ©  and 
its  tail),  41.5%  of  the  NE  capacity  must  have  been  employed  in  the 
discharge  regime.  On  the  other  hand,  due  to  the  irreversible  loss  of 
Li  inventory  typically  occurs  in  the  first  cycle  of  the  formation,  most 
of  the  cell  design  is  loaded  with  NE  in  excess  of  capacity  in  order  to 
capture  the  full  capacity  of  the  PE.  As  such,  the  PE  loading  usually 
is  the  one  limiting  the  cell  capacity.  Therefore,  after  the  initial  SEI 
formation,  only  91 .5%  of  the  capacity  (41 .5%  from  O  and  nearly  50% 
from  the  remaining)  in  the  graphite  is  utilized.  We  thus  conjecture 
that  the  loading  ratio  of  NE/PE  is  around  1.08. 

Fig.  9(a)  and  (b)  presents  the  IC  curves  associated  with  different 
cycling  rates.  At  the  first  glance,  they  all  appear  similarly  in  behav¬ 
ior,  except  for  5C,  which  exhibits  a  significant  change  in  behavior 
during  discharge  and  almost  no  signature  during  charge,  likely  due 
to  a  high  degree  of  polarization.  In  order  to  compare  IC  curves  of 
different  rates,  it  is  easier  to  exhibit  them  in  terms  of  IR-free  voltage 
instead  of  the  measured  cell  voltage.  The  polarization  resistance  in 
the  charge  and  discharge  regimes  is  thus  needed.  Fig.  10  presents 
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the  evolution  of  the  cell  voltage  as  a  function  of  the  IR  drop  at  (a)  the 
BOD  and  (b)  the  beginning-of-charge  (BOC)  at  different  rates.  Since 
all  rates  have  almost  the  same  starting  SOC,  the  initial  IR  drop  can  be 
used  to  determine  the  cell  internal  resistance  at  both  BOD  and  BOC 
Only  rates  lower  than  1C  were  used  to  ensure  that  the  observation 
is  still  in  the  ohmic  regime  of  the  Tafel  behavior,  where  the  voltage 
response  to  current  is  linear.  The  polarization  resistance  was  found 
to  be  132  m^  Ah  and  216  m^  Ah,  or  69  m^  and  114m£2,  respec¬ 
tively,  at  BOD  and  BOC  for  the  NSC.  It  indicates  that  the  resistance 
at  the  BOC  is  significantly  higher  than  that  at  the  BOD. 

Fig.  9(c)  and  (d)  presents  the  IC  curves  without  polarization- 
induced  voltage  drop  for  a  better  comparison  of  the  IC  peaks  in  the 
discussion.  The  rate  effect  on  the  cell  behavior  can  be  separated 
into  two  parts:  one  for  the  rates  below  1  C  and  the  other  above.  For 
the  lower  rates  only  peak  broadening  is  observed,  indicating  that 
the  solid-state  diffusion  remains  as  the  rate-limiting  factor.  The 
peaks  in  the  low  voltage  region  of  the  discharge  regime  and  the 
high  voltage  region  of  the  charge  regime  appear  more  broadened 
than  those  in  the  other  regions,  likely  a  result  of  the  grain  size  (long 
diffusion  path). 

Substantially  different  behavior  in  the  IC  peaks  was  observed 
for  the  higher  rates.  It  seems  that  all  IC  peaks  with  O  are  fading 
in  rate  capability  as  the  rate  increases.  This  phenomenon  is  not 
peculiar,  as  it  is  known  that  reaction  O ,  which  corresponds  to  the 
transformation  of  LiC6  to  UC12,  is  kinetically  slower  than  other 


anodic  reactions  leading  to  the  lithium  depletion  in  the  graphite 
structure.  Therefore,  the  peaks  associated  with  O  (thus,  O  *(©  , 
©  ,  and  ©  ,  respectively))  are  broadened  more  than  the  others  as 
the  rate  increases  and  start  to  overlap  with  peaks  associated  with 
reaction  ©  ,  especially  ©  *©  .  This  behavior  leads  to  a  longer  time- 
resolved  extent  of  reaction,  as  on  the  surface  of  the  active  grains, 
the  cell  voltage  indicates  that  reactions  ©  *((D ,  ©  respectively)  have 
already  occurred,  whereas  in  the  core  the  composition  remains  at 
O  *(©  and  ©  ,  respectively)  and  needs  to  catch  up  with  the  sur¬ 
face  reaction.  An  extreme  case  is  observed  in  5C,  where  reaction 
O  *©  occurs  slowly  throughout  the  grain,  while  the  cell  voltage 
as  measured  on  the  grain  surface  is  approaching  reaction  ©  .  This 
phenomena  and  their  consequence,  as  exhibiting  on  IC  signals, 
continue  to  play  out  in  the  cell  reactions  in  a  dynamic  manner, 
explaining  the  extremely  broadened  peak  in  convolution  at  this 
rate. 

As  a  general  trend,  it  appears  that  LiMn^Ni^Co^C^  is 
affected  by  the  increase  of  current  (rate)  more  than  LiMn204  (thus, 
worse  rate  capability).  The  evolution  of  peaks  O  *©  and  O  *© 
explains.  The  transition  from  ©  to  ©  occurs  on  the  staging  plateau  O 
,  which  eliminates  the  contribution  from  the  NE  to  the  change  of  the 
peak  shape  over  cycling.  Upon  cycling,  peak  O  *©  is  clearly  broad¬ 
ened,  delaying  the  onset  of  the  next  reaction;  whereas,  the  shape  of 
peak  O  *©  appears  only  slightly  affected.  This  observation  affirms 
that  LiMn204  is  a  constituent  that  has  higher  rate  capability  for 
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has  profound  potential  to  analyzing  the  degradation  of  the  cell  in 
the  second  part  of  this  study.  This  implication  is  important  to  us, 
as  these  cells  are  intended  to  support  the  design  of  an  experiment 
that  uses  a  large  matrix  of  cells  and  test  conditions  to  study  path 
dependence  of  cells’  degradation  under  PHEV  duty  cycles  and  con¬ 
ditions.  Between  UH  and  INL  more  than  60  cells  are  under  testing 
with  different  aging  protocols.  The  result  of  this  study  on  cell  varia¬ 
tions  established  a  reliable  basis  for  comparison  of  future  aging  test 
results.  It  also  helps  to  establish  a  consistent  diagnostic  methodol¬ 
ogy  to  enable  a  higher  possibility  of  success  in  this  experiment  and 
to  assure  that  the  observed  discrepancies  in  the  test  data  can  be 
attributed  to  aging  protocols  with  certain  confidence. 
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Fig.  10.  Evolution  of  the  cell  voltage  due  to  IR  drop  as  a  function  of  rate  at  (a)  BOD 
and  (b)  BOC  in  cell  #1. 


high  power  delivery,  whereas  the  LiMrq^Ni^Ccq^C^  contributes 
more  to  the  capacity. 

Change  in  behavior  of  the  IC  peaks  over  aging  shall  give  us  a 
vantage  point  to  infer  mechanisms  that  affect  capacity.  Seeing  how 
active  materials  and  their  constituents  change  over  time  and  aging 
conditions  will  be  informative  toward  improving  cell  design  and 
choice  of  materials,  as  well  as  exploiting  the  best  conditions  for 
which  to  use  these  materials.  This  information  will  be  incorporated 
into  physics-based  degradation  models  to  elucidate  further  knowl¬ 
edge  for  root-causes  of  cell  aging  and  will  be  discussed  in  the  part 
II  of  this  work. 

5.  Conclusion 

A  batch  of  ten  1.9  Ah  commercial  18650  cells  was  investigated 
to  characterize  its  baseline  performance  and  intrinsic  cell-to-cell 
variations.  It  was  found  that  the  cells  are  of  high  quality  with  little 
cell-to-cell  variations,  thus  establishing  a  high  confidence  level  for 
analyzing  and  interpreting  test  results  from  a  large  matrix  of  cells 
that  undergo  various  test  conditions.  Moreover,  these  cell  varia¬ 
tions  should  be  incorporated  and  accommodated  into  future  cell 
and  pack  performance  models  and  simulations  to  enhance  accuracy 
in  prediction. 

The  analysis  of  the  incremental  capacity  (IC)  curves  showed  that 
the  complex  IC  signature  could  be  deciphered  and  that  information 
on  each  specific  active  component  in  the  cell,  LiMn1/3Ni1/3Co1/302 
and  LiMn204,  in  the  PE  and  graphite  in  the  NE,  could  be  gathered 
respectively.  From  the  analysis  of  the  relative  IC  peak  capac¬ 
ity  in  some  of  the  reactions  specific  in  each  components,  we 
are  able  to  conjecture  that  the  c-PE  comprises  about  2/3  of 
LiMn1/3Ni1/3Co1/302  and  1/3  LiMn204,  and  similarly  that,  after  the 
initial  SEI  formation,  the  loading  ratio  of  NE/PE  is  about  1.08.  This 
inference  ability  to  deciphering  IC  behavior  in  the  cell  chemistry 
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